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ESTIMATION OF THE LENGTH ANI) ORIENTATION OF THE LINE 
BETWEEN TWO CLOSELY CO-ORBITING SATELLITES 


John Urn star* 
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Abstract 


The problem of estimating the baseline length and 
orientation using a combination of on-board sensor data 
and processed ground tracking data is considered. A 
minimum number (3) of simple on-board sensors are 
assumed. Some of the variables are assumed to be esti- 
mated from ground tracking information only. These 
estimates become part of the observations, along with 
the on-board sensor outputs, for estimation of the re- 
maining variables. A linear estimate of the remaining 
variables is derived via a Kalman filter. This sequen- 
tial processing obviates the need for the simultaneous 
processing of ground track and on-board data. The im- 
portance of deciding which variables are to be estimated 
from the ground and which are to be estimated using the 
combined data is discussed. An example demonstrates 
that baseline length and orientation can be estimated 
wit v in a few feet and fractions of a degree. 

INTRODUCTION 

Recently interest has been expressed in using pairs 
and even clusters of satellites in close proximity to per- 
form a variety of space missions. 1,2 One particular 
area in which dual satellite systems should prove very 
useful is in low frequency radio astronomy. A pair of 
satellites may be used as an interferometer to synthesize 
a large aperture radio telescope.- 4 * 0 As the pair moves 
in orbit the location and frequency characteristics of 
various radio sources may be determined. 

The satellites in these systems may be tethered 
thereby forming a very large "dumbbell” satellite, or 
they may be co-orbiting** but physically separate. In 
either case, however, when the pair of satellites is used 
as an interferometer, the length and orientation of the 
line between them (the baseline) must be determined to a 
relatively high degree of accuracy. 

The problem considered here is the estimation of the 
baseline length and orientation of two separate but co- 
orbiting satellites. The basic assumption is that the 
baseline is small compared to the orbit dimension. For 
example, the radio astronomy mission may call for an 
altitude of 8000 mJes and baseline of only one or two 
miles. Because of the closeness of the two satellites the 
baseline length and orientation cannot be adequately esti- 
mated from ground based measurements alone. Small 
errors in estimates of the position vectors of each satel- 
lite reflect as large errors in the baseline orientation 
estimate. Thus there is a requirement for some on-board 


*Aerospace Engineer 

** Co-orbiting is taken to mean the two satellites are in 
the same nominal orbit and the difference in the times 
from perigee is a constant. 


sensors. The baseline length and orientation could be de- 
termined with an intersatellite ranging sensor and an 
optical sensor on one satellite viewing the second one 
against the star field background (the latter sensor de- 
termining baseline orientation). In this case on-board 
sensor data alone would solve the problem. An optical 
sensor of this type, however docs not presently exist. 

The approach here was to use simpler, state-of-the-art 
sensors. These sensors however do not provide enough 
information to make the system observable. The solution 
of the problem therefore lies in the combination of the 
on-board sensor data and ground tracking data* 

The problem is formulated in terms of the small 
perturbations from the nominal circular orbit common to 
both satellites. The orientation of the baseline may be 
described by two angles, one in the nominal orbit plane 
and one between the baseline and this plane. The linear- 
ized equations for the in-plane and out-of-plane motions 
are not coupled. The problem is therefore reduced to 
two simpler ones. The state variables are defined as the 
sums and differences of the satellite positions and veloc- 
ities. Some of the variables are estimated using processed 
ground based measurements. These estimates, along with 
the on-board sensor information, are used to estimate the 
remaining variables. The estimates of the two angles and 
baseline length are derived from the state variable 
estimates. 


PROBLEM FORMUIATION 

Assuming a circular earth the inertial frame to be 
used as a reference is the XYZ frame where X and Y 
define the nominal orbit plane and the origin is at the 
center >f the earth (Figure 1). 

The orientation of the baseline relative to this frame 
may be described by the angles ^ and /. Beta is the 


t 



Figure 1. Orbital Configuration of the 
Dual Satellite System 






mul<* bHwivn t hi* Baseline ami Uu* X-Y plane ami . in 
the- angi** Between the Baseline planar projection and the 
X ums (figure h, These angles and tin* baseline length 
h ma\ oe determined Irom Uu* variables shown in the 
figure, rims the estimates ot , and B are functions 
el the vectors locating the satellites A and B In the* XYZ 
Irnme. 


2, A \aw sensor measuring the* A satellite yiw angle 
(the attitude angle around the o 1 it radius 
\ ector), 

a. An optica) "H-tracker" locating the B satellite 
with respect to the A satellite altitude* reterenee 
frame 1 . This measu res the angle , 


The equations of motion of each satellite with respect 
to this frame aria 


The latter two angles are shown in Figure* 2, Both , , 
arc 1 nominally zero. The sensor is assumed to Ik* in- 
despondent of tin 1 small variation in . 
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In equations la, I), e cylindrical coordinates are used 
where U and arc 1 in the X-Y plane and z is the dis- 
tance out of the plane. The gravitational constant is 
The state variables for the A satellite* are defined as 
components of the* vector q. 


Figure 2. Definition of , , , 
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The relation between the state variables and the ob- 
servables C y ( i ) ) may be determined from the geometry 
of the problem. As a function of the state variables, I) is 


The state variables for satellite* B are similarly defined 
as components of the vector p. Components of the vectors 
q and p are indicated by numerical subscripts. The nom- 
inal orbit radius is R 0 . When the equations are linearized 
with respect to the nominal circular orbit the system 
equations are: 


b U\f , p* - 2q 1 p | cos 
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The first variation in b is therefore 
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In equations (2) and (3) *q(t) and ip(t) are the 
first variations of q(t) and p(t). The dynamic noise 
is u(t) (the same for lx>th satellites because of their 
proximity) and the sensor noise is v(t). The state 
transition matrix <I>(t,t 0 ) is listed in the appendix. The 
lack of coupling between the planar variables, -q, , . p t 
through oq 4 , p 4 and the non-planar variables oq s , op 5 , 
,q 6> p 6 is evident. 

To determine the relationship between the output 
vector cy(t) and the state vectors the matrix H(t) must 
be defined. This requires definition of the sensors. The 
following sensors are assumed to be located on A. 


b Cdi < Pi) sir M0rh (‘'M 2 'P 2 ) cos 


where 6 c is the angle subtended by the nominal baseline 
B 0 when A and B are in the nominal oi’bit, i.e.,sin( 2) 
- B 0 /2R q . Equation (5) is independent of the z variables 
q s , p s . The intersatellite range sensor is the only 
on-board sensor required for the estimation of b and . 
That the estimate of b may be derived from this infor- 
mation is obvious. It will be shown that the > estimate 
may also be derived using only the intersatellite range 
information (b measurement) and processed ground 
tracking information. 

The angle p is expressed as a function of the state 
variables as follows 


(p 5 q 5 ) 

/n 2 - p 2 \ 

qf-P?-2q lPl co»^_j,(p s q 5 ) 2 


This angle is not directly observable but is related to the 
angles p (Figure 2): 


1. An inter sate llite range sensor measuring base- 

line length b. 


♦Brackets [ ] indicate a matrix 
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(7) 


Planar Solution 


Using equations (<») and (7) the iirst variations in the ob- 
servables . and - are i> t , and 

"r • ( P 5 q S ) ( H ) 

u o 

The introduction of requires the introduction of the 
linearized yav attitude dynamics 

oH»(t) ^(t.t 0 ) b»(t 0 ) * w(tj v 9) 

where o V ( t ) is the rxl state vector with ■ v as its first 
component and w(t) is the dynamic noise. The corre- 
sponding transition matrix is <l> 3 (t, t 0 ). 

Examination of equations (5) and (8) shows that planar 
and non-planar observations are uncoupled just as the dy- 
namics. The problem may therefore be considered as 
two uncoupled problems. Equations (2) and (3) hold for 
the planar problem (but only the first four components of 
oq and .;p are used andOjft, t 0 ) replaces 4>(t,t 0 )) In 
the planar case .y(t) = it b( t ) , 11= fsin(A^/2), 
cos(A;y2), 0, 0, sin(Atf/2), -cos(Afc/2), 0, 0] and v t (t) 
is the range sensor noise. The tilde (~) refers to the 
observable with measurement noise included. 

The corresponding equations for the non-planar 
problem are equations (10), (11), and (12). 


Estimates of the variables b ami • may be determined 
as functions of the estimates of the planar state variables. 
These variables must therefore be related to the state vari- 
ables. For b the relationship has already been shown to be 
equation (5). The first variation in ■ is 

f-n, p, ) . -ji- * p 2 ) (i») 

As the problem stands there are eight variables 
. p t , i-l,4 and one observable cb. With the single on- 
board sensor the system is not observable. Not enough 
information is available to allow an estimate of and 
p. This problem may be overcome by redefining the 
state variables as sums and differences of the cq t ’s and 
op^s and then estimating some of the new variables from 
ground observations. 

Let 

f oq k - Pj i 1, 6 (14a) 

' < 

hq u6 1 °Pl-6 1 7 - 12 ( 14l >) 

With this change of variables the planar system equations 


become 
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In order to provide the necessary information to 
make the system observable it is now assumed x 2 is 
estimated from ground observations (Figure 3). Thus 
equation (16) is no longer needed. Note also that ib con- 
tains components from both x 1 and x 2 . 

(11) 

Denote the estimate of x. by x ; and the error in the 
estimate by e r 

Xj = X. _ e . i r 1, 12 (18) 

Substituting for x 7 in equation (17) results in 


H, = 


0 

0 


r 


0, I - 1/Bg 0 



The and rj sensors have noise components v 2 (t) 
and v 3 (t) respectively. The components of v(t) a re 
all assumed to have zero mean. 


I>b -XjCOS (\~y) + X 7 sin (^r)~ e 7 Sin (^) + V l (t) (19) 

At any time the estimate X ? is a constant entered from 
an external source (the output of the x 2 estimation 
process). Thus it may be incorporated in the measure- 
ment &b. The statistical properties of e 7 and v t (t) 


♦The T indicates transpose. 
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PROCESSING OF 
GROUND TRACK 
DATA 


*2 


PROCESSING OF 
COMBINED 
DATA 


9 , 0 , b 


In equation (2.1) E * * is the expected value operator 

and • * , i arc the variances o! the range sensor noise 
and x, estimate error respectively. Since v, ft) and 
are uncorrelated E (v, ft ) . - ' 0, 


B> using equations (l-"o and (14) and x l and c r 
may In* estimated 


1 . 
— x. 


2R„ 


(2«i) 


In equation (20) x, is from x l , x s is from x 2 and * . is 
the estimate of^the variation b > . The error of the esti - 
mate is r y The variance of the estimate 

error is 


Figure 3. Information Flow Diagram 


are assumed to he known and e ? and v t ft) are uncor- 
related. Taken together these terms constitute the total 
measurement noise. The output equation therefore re- 
duces to 


■•1. - 


/ ' 

\ 

0. 4-os (_ 

) , 0. 0 

\ 
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sin 



Vj ft ) (20) 


Equations (15) and (20) are in the correct form to allow 
x 1 to he estimated. 

The vector x 1 may be estimated by means of the 
Kalman filter equations • 

x (t + 1') - 'Iq (t + 1, t ) x ft ) + K (t + 1) (b y (t + 1) 

- lift + 1) <;>! (t 4 1, t) X(t)) (21) 


Kft + li = i\ ft + l) n T ft + n ■ 

( H ( t 4 1) p t (t 4 l)H T (t 4 1) + R(t 4 l))- 1 (22) 


,2 ,2 
• 2 =_i+ JL 


R 0 B o 


(27) 


where v 2 , c* are tiie variances of the x, and x 8 esti- 
mate errors and >• is the correlation coefficient between 
these estimate errors. Because of the separate process- 
ing of x 1 and x 2 , ,, is not directly available. However 
!’ 2 is hounded by 

y * 


a 


2 

y * 


!l + iL , tr > h 

B 2 4 R 2 R o B o 


(28) 


Because of the basic assumption that B 0 « R 0 it will turn 
out that „ = < r y . This will be demonstrated in the 
example. 


The estimate of eh and variance of the estimate 
error arc 


eb = H x 1 = cos 



(29a) 


2 _ 
b ~ 


2 

C O S ^ 



(29b) 


The covariance matrix of the estimate error P (t 4 1) is 
determined by equations (23) and (24). 

P t (t + 1) = <\\ ft 4 1, t) P(t) (t 4 1) 4 Q(t) P 3 ) 

P(t 4 1) =P t (t 4 1) - K (t + 1) H T (t 4 1) P t (t 4 1) (24) 

In equations (21) through (24) x (t + l)is the estimate 
given measurements through time t + 1, x (t) is the esti- 
mate given measurements through time t . The measure- 
ment interval has been normalized to one in these equa- 
tions. The covariance matrices of the dynamic noise and 
measurement noise are Q(t) and R(t) respectively. 

For the planar problem 



Non-Flanar Solution 

The angle /i is defined as a function of the state 
variables in equation (6). With the variables defined in 
equation (15) the first variation is 

bfi = - i- X . (30) 

B o 3 

Thu i the estimate of fi is a function of the non-planar 
variable x $ only. 

The system equations for this problem (using equa- 
tions (10), (11) and (12)) are 
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Once an estimate of x s is available (via the filter equa- 
tions) the estimate of ifi and variance of the estimate 
error are respectively 
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EXAMPLE 


(33a) 


(33b) 



Consider a dual satellite system with the mission of 
low frequency radio astronomy. The co- orbiting satel- 
lites are the receivers of an orbiting interferometer. In 
order to determine the location and frequency character- 
istics of the emitting sources the baseline length (b) and 
orientation (y, /?) must be known. 

A digital simulation was carried out using the follow- 
ing parameters 

R 0 = 70 • 10 6 feet 

a) Q s 2.03 > 10" 4 sec -1 

B 0 = 10 4 feet 

r = 100 sec (sampling interval) 

The dynamic noise was assumed to be zero. The meas- 
urement noise was assumed to be white Gaussian with 
zero mean. 

cr vl = 5 feet 

? v2 -cr v3 =8.7 milliradians (0.5°) 



Figure 5. Estimate Error and Standard Deviation of x 2 


The estimate errors of x 7 and x g after processing of 
ground based measurements were assumed to have 
standard deviations 

cr 7 = 1500 feet 

o-g = 300 feet 

Sample runs for the planar case are shown in Fig- 
ures 4 and 5. These figures show the estimate errors of 
x ( and x 2 (i.e. e, and e 2 ) and their respective standard 
deviations. Good estimates of Xj and x 2 are available 
after 50 samples (5000 second?). Using o-j and <r 2 after 
50 samples, 35 feet and 2.5 feet respectively, in equations 
(28) and (29b) results in cr y , = 35 * 1.0 ~ A rad. (0.2°) and 
cr b = 2 fi feet. Because or the ^sumption B 0 «x R 0 the 


second and third terms of equation. (23) are negligible. 
Thus Vy. - <r y - CT t /B 0 . For the non-planar case the 
sample run is shown in Figure 6. A simple second order 
system was assumed for the yaw attitude dynamic s. The 
resulting a after 50 samples is 25 • 10” 4 radians (0.14°). 

DISCUSSION 

Two points should be emphasized with respect to the 
preceding results. 

There may be cases in which the B estimate is not 
required either because the out-of-plane motion is not 
significant or because the sources cf interest are in or 
near the orbit plane. In this case the quantities of 
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The* three parameters describing the length and or i 
ontatlon of the line lietweon two closely co-orbiting nat«‘l 
lites cannot he adequately estimated from ground mean 
uromonts alone because small uncertainties in tlu sate lilt 
location vectors contrilxitc to large uncertainties in the 
l>aseline orientation. A minimum numlx*r of simple on- 
board sensors alone cannot do the job completely Ijeoauso 
they do not provide enough information in themselves, 
However, by combining the information provided by three 
on-board sensors (two angle sensors and an Intersatellite 
ranging system) with the processed ground tracking in- 
formation, a good linear estimate of these three param- 
eters can be obtained. Care must Ik* taken in defining 
the variables to Ik* estimated to insure a low error in the 
final output. 


Figure <i. Estimate Krror and standard Deviation ofx s 
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APPENDIX 


The state transition matrix <1> (t, t 0 ) is 
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